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ABSTRACT 

The giant radio galaxy M87 with its proximity (16Mpc), famous jet, and very massive black hole ((3 - 
6) X 10 J M o ) provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) 
y-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has 
been established as a VHE y-ray emitter since 2006. The VHE y-ray emission displays strong variability on 
timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the 
MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected 
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triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent 
sampling of the VHE y-ray light curve enables one to derive a precise temporal characterization of the flare: the 
single, isolated flare is well described by a two-sided exponential function with significantly different flux rise 
and decay times of r™ e = (1.69 + 0.30) days and r | LCay = (0.611 ± 0.080) days, respectively. While the overall 
variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 
2008, they share very similar timescales (~day), peak fluxes (flCojsTeV — (1 - 3) x 1CL 11 phcm _2 s _1 ), and 
VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in 
contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided 
with a VHE flare. On the other hand, Chandra X-ray observations taken ~ 3 days after the peak of the VHE 
y-ray emission reveal an enhanced flux from the core (flux increased by factor ~ 2; variability timescale < 2 
days). The long-term (2001-2010) multi-wavelength (MWL) light curve of M 87, spanning from radio to VHE 
and including data from 1 1ST, LT, VLA and EVN, is used to further investigate the origin of the VHE y-ray 
emission. No unique, common MWL signature of the three VHE flares has been identified. In the outer kpc 
jet region, in particular in HST-1, no enhanced MWL activity was detected in 2008 and 2010, disfavoring it 
as the origin of the VHE flares during these years. Shortly after two of the three flares (2008 and 2010), the 
X-ray core was observed to be at a higher flux level than its characteristic range (determined from more than 
60 monitoring observations: 2002-2009). In 2005, the strong flux dominance of HST-1 could have suppressed 
the detection of such a feature. Published models for VHE y-ray emission from M 87 are reviewed in the light 
of the new data. 

Subject headings: galaxies: active - galaxies: individual (M 87) - gamma rays: galaxies - galaxies :jets; nuclei 
- radiation mechanisms: non-thermal 
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1. INTRODUCTION 

The giant radio galaxy M 87 provides a unique environment 
to study relativistic plasma outflows and the surroundings of 


super-massive black holes (SMBH). Its prominent jet (Cur 
tis]|1918) is resolved from radio to X-r ays displaying com- 


plex st ructures (knots, diffuse emission; Perlman et al.|1999 


2001a), strong variability ( |Harris et al.[ 2003 20 06[>, and ap 


parent super-luminal motion ( Biretta et al.|1999[[Cheung et al. 
2007). With its proximity (16.7 + 0.2 Mpc; |Mei et al.|2(. 
and its ver y mass ive black hole of A7p, n - (3 - 6) x 1 0 9 M 
(Macchetto et al. |1997 |Gebhardt & Thomas|2009[ ) high res- 
olution very long baseline interferometry (VLBI) a 


very long baseline interferometry (VLBI) at radio 
wavelengths enables one to directly probe structures with 
sizes down to < 200 Schwarzschild radii. From the detection 
of super-luminal features in the jet in the optical and radio the 
jet orientation angle towards the line of sight at the sub-kpc 
scale is limited to 6 < 20° (Biretta et al.| 1999) Cheung et al. 
[2007) . 

Evidence for very high energy (VHE; E > 100 GeV) y-ray 
emission from M87 was reported by the HEGRA collabo- 

and was later con- 


ration in 2003 (Aharonian et al. 2003 
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firmed by H.E.S.S., VERITAS, and MAGI C (Aharonian et al. 
|2006a[ |Acciari et al.||2008) | Albert et al.||2008a| ). While the 
large majority of active galactic nuclei (AGN) with detected 
VHE y-ray emission are strongly beamed sources, M 87 is one 
of only four known VHE AGN with weak or, at most, mod- 
erate bea ming; the other three being the radio galaxies Cen - 
taurus A (Aharonian et al.|2009), IC 310 (| Aleksic, et al.|2010) , 
and NGC 1275 (Mariotti & MAGIC Collaboration|2010a) . In- 
terestingly, for such a weakly beamed source, M87 shows 
strong variabi lity at VHE with timescales of the order of d ays 
(Aharonian et al.||2006a[ |Albert et al.||2008al |Acciari et al.| 
2009). This points, through the causality argument, towards a 
compact emission region < 5 X 10 1S (5 cm (5 being the Doppler 
factor of the emitting plasma) corresponding to only a few 
Schwarzschild radii Rs = 2 GMbh/c 2 - 10 15 cm. At GeV en- 
ergies M 87 has recently been detected as a weak source by 
Fermi-LAT (Abdo et al.|2009) . 

The exact location of the VHE y-ray emitting region in 
M 87 remains elusive. The angular resolution of ground-based 
VHE instruments is of the order of 0.1° (corresponding to 
~ 30kpc projected size) and, therefore, does not allow for a 
direct precise determination of the VHE y-ray emission site in 
the inner kpc-scale structures, although the outer radio lobes 
can be excluded as the origin (Aharonian et al. 2006a). To 
further investigate the location of the VHE y-ray emission site 
and the associated production mechanisms, variability studies 
and the search for correlations with other wavelengths have 
successfully been utilized (e.g. Acciari et al.~| [2009 ). Of par- 
ticular interest are radio observations, since they allow for 
the highest angular resolution, and X-ray observations, due 
to their potential connection with the VHE y-ray emission in 
e.g. synchrotron self-Compton (SSC) models. 

Up to now, three episodes of enhanced VHE y-ray emis- 
sion have been detected from M 87, with details on the lat- 
est one, observed in 2010, being reported in this paper. The 
first one, detected in 2005 (Aharonian et al.|2006a) , coincided 
with a n extreme multi-frequency ou tburst of the jet feature 
HST-1 (Harris et al.|2003||2006) which has also been dis- 


In the fo llowing Mrh - 3 x IQ 9 M© is adopted 


cussed as a possible VHE y-ray emission site (e.g. Stawarz 
et al.||2006 Cheung et al.||2007 Harris et al. |2009 1 . During 


the second flaring episode, detected in 2008, HST- 1 was in 
a low flux state, but radio measurements at 43 GHz with the 
Very Long Baseline Array (VLBA) showed a flux increase in 
the core region within a few hundred Schwarzschild radii of 
the SMBH, suggesting the direct vicinity of the SMBH as the 
origin of the VHE y-ray emission (Acciari et a l.|2009 1 . This 
conclusion was further supported by the detection of an en- 
hanced X-ray flux from the core region by Chandra. 

In this paper results from a joint observation campaign of 
M 87 in 2010 are presented and discussed in the broader con- 
text of the multi-wavelength (MWL) behavior of M 87 over 
the past ten years. During the campaign, a high flux state at 
VHE was detected ( Mario tti & MAG IC Collaboration|2010bl 
[Ong & Mariotti ||20 10) . Characteristics of the VHE flare are 
investigated and possible correlations with other wavelengths 
are discussed. New observational results from H.E.S.S., 
MAGIC, VERITAS, Fermi-LAT. Chandra , HST, LT, VLBA, 
MOJAVE, VLA, and the EVN are present ed | 1U7 | 

In Sec. [2] the instruments and the data are introduced. In 
Sec. [3] the characteristics of the VHE high state are investi- 
gatedand compared to previous flares. New results from opti- 
cal polarimetry observations with the Hubble Space Telescope 
and the Liverpool Telescope are presented and discussed in 
Sec. [T] In Sec. [5] the VHE results are confronted with the 
broader MWL picture and the theoretical implications con- 
cerning the VHE y-ray emission site are presented in Sec. [6] 
The paper concludes with Sec. [7] 

2. DATA 

2.1. Very high energy (VHE) 

M 87 has been the target of several coordinated VHE mon- 
itoring campaigns since 2008. In 2010 the source was jointly 
monitored by MAGIC and VERITAS. During the campaign 
[Mariotti & MAGIC Collaboration) (2010b| reported an in- 
creased flux from the source in February 2010, but follow-up 
observations did not reveal further activity at VHE. A sec- 
ond increase in flux, marking the onset of a strong VHE flare, 
was detected in April 2010 by MAGIC and VERITAS ( Ong &| 
|Mariotti |2 0 10) triggering further ToO (Target of Opportunity) 
observations by H.E.S.S., Chandra, the VLBA, and the EVN. 

For the 2010 campaign, M87 has been observed for a total 
of ~80h from December 2009 to June 2010. The data from 
different VHE instruments have been combined after separate 
analysis within the individual collaborations. Integral fluxes 
for the VHE band <J> V hk are calculated above an energy of 
350 GeV. Observations taken or published with a different en- 
ergy threshold are extrapolated using the average measured 
VHE spectrum, which is well described by a power law spec- 
trum dN/dE oc E 1 with photon index T = 2.3 ±0.11 ( Albert) 
|et al.|2008a) . 

Given the indications for spectral variability of M87 at 


spi _ 

VHE (Aharonian et al.|2006a[ [Albert et al.|2008a) |Aliu et aL 


105 See also 


Le Bohec et al. 


(2004 


The outburst was also followed by the ejection of apparent superlumi- 


201 1 ) using a single photon index for the flux extrapolation 
could introduce a bias in the light curves. For the flaring 
states the harde st sp ectral index reported is T = 2.21 +0.18 
(Albert et al. 2008a) and, therefore, the systematic error in- 
troduced when using T = 2.3 is completely negligible com- 
pared to the typical statistical error of 10-20% of the flux dur- 
ing such states. For the quiescent state spectral indices up 

nal radio components from HST- 1 (Cheung et al.|2007|. 

107 Full names of the instruments can be found in Sec.|2] 


2010 VHE flare & 10 years of MWL observations of M 87 


5 








3.5 
3 

2.5 
2 

1.5 
1 

0.5 

0 

-0.5 

7 5 

7 e 4 
“ 3 
S 2 
£ 1 
0 

1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

. 12 
' 10 
! 8 
i 6 

5 4 

' 2 
0 

0.8 

0.6 

0.4 

0.2 

0.0 

2.0 

1.5 


0.5 

0.0 

3.5 

3.0 

2.5 

2.0 

1.5 
1.0 
40 
35 
30 
25 
20 
15 
10 
05 
00 


Year 

^ ^ T?* TjfP ^ 



52000 


52500 


53000 


53500 54000 54500 55000 55500 

MJD 


Fig. 1. — Multi-wavelength light curve of M 87 from 2001 to 201 1. The VHE y-ray flux (top panel) is calculated above an energy threshold of 350 GeV (see 
text). Separate fluxes for the core and HST-1 are shown in cases where the instrument resolution is sufficient to separate the two components. Gray vertical bands 
mark the times of increased VHE activity in 2005, 2008, and 2010 (see Fig. IT). The dashed line and the gray horizontal band in the 2nd panel marks the average 
flux with 1 s.d. error measured by Fermi-hAT . The radio flux of HST-1 at aifferent frequencies has been normalized to the 5 GHz flux assuming a spectrum 
S v ~ v~ a with a = 0.6. All flux errors shown are the 1 s.d. statistical errors except for the LT data where the uncertainty on the contribution from the galaxy is 
included in the error bars. For details on the data, the data analysis, and references see text. 


to r = 2.60 ± 0.30 have been reported ( Albert et al.|2008a) . 


This could, in principle, create a bias for the H.E.S.S. data set, 
where the flux is extrapolated down to 350 GeV from energies 
above 500 GeV. For example, for spectral indices of Tj = 2.3 
and [A = 2.6 the flux extrapolated down to 350 GeV from 
700 GeV would differ by a factor (350 GeV/700 GeV)' 2-1 ' ~ 
0.8, which would imply errors of order ~20%. Given that 
the typical flux error for the nightly averaged flux bins in the 


Given recent results on the Crab Nebula (the reference 
source for ground-based VHE instruments) indicating only 
small systematic offset between t he en ergy scale of differ- 


ent instruments ( Meyer et al.||201 0p^ 


and the general good 
agreement of the flux measurements during quasi simultane- 
ous observations, the systematic error between the different 
instruments is estimated to be small compared to the statisti- 
cal error of individual measurements of the 2010 campaign. 


quiescent state is of order ~ 1 00 % (i.e. the individual nightly 

Additional archival data from 

Aharonian et al. (2006a), 

Ac- 

flux points are not significant detections) the systematic un- 

ciari et al. (12008), Albert et al. 

2008a), Acciari et al. (2010), 


certainty introduced by using a single photon index for ex- 
trapolation can also safely be neglected. 


108 Energy scale scaling factors relative to the Fermi - LAT energy scale of 
0.961 ± 0.004 for H.E.S.S. and 1.03 ± 0.01 for MAGIC are derived. 
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and |Aleksic et al.| ((201 laf are also shown in the light curves 
(Fig- [l] and Fig. |2]». 

In the following, the main characteristics of the VHE ob- 
servatories involved in the M 87 campaign are reviewed, and 
details on the corresponding data-sets are presented. 


H.E.S.S. — The 
(H.E.S.S. 


109 


High Energy Stereoscopic System 
is a system of four large (13 m mirror di- 
ameter) imaging atmospheric Cherenkov telescopes (IACTs) 
for the detection of VHE y-rays, located in the southern 
hemisphere in Namibia (23°16'S, 16°30'W; 1800m above 
sea level). It has been in operation since 2002, with the full 
array completed in 2004 ( |Hinton||2004) >. H.E.S.S. observed 
M87 for 10.7 h in 2010 (dead-time corrected) yielding a 
total detection signific ance of 9.7 standard deviations (s.d.; 
following |Li & Ma| 1983[> using a standard H illas-type 
analysis and cuts fromjAharonian et al.||2006b| (software 
version hap- 18-06). M 87 culminates at ~35° zenith angle 
at the H.E.S.S. site resulting in an energy threshold of 
Ethr > 500 GeV for the analysis and data-set. The integral 
flux is extrapolated down to 350 GeV using the average 
energy spectrum (see above). In addition, results from a 
re-analysis of the 2004 and 2006 H.E.S.S. data in nightly 
flux bins, utilizing the same analysis as discussed above, 
are presented in the light curves. Cross check on the results 
with data from an independent calibration chain and utilizing 
different analysis have been performed and good agreement 
is found. 

MAGIC — The M ajor Atmospheric Gamma -ray Imaging 


Cherenkov (MAGICl 110 telescope system consists of two 
17 m diameter IACTs located at the Roque de los Muchachos 
Observatory, on the Canary Island of La Palma (28°46' N, 
17°53' W; 2200 m above sea level). Since 2005 M 87 has been 
regularly observed by MAGIC with a single telescope (|Albert| 
|et al.|[2008b| jAleksic et al.||20lTa) . In the Autumn of 2009 
the system became stereoscopic, with the commissioning of 
a second telescope and the stereo trigger, resulting in an al- 


most doublin g of its sensitivity ( Cohn et al. |2009| |Aleksic 


|et al.||201 lb[ ). In 2010, M87 observations were conducted, 
for the first time, in stereoscopic mode. 20 h of good quality 
data were taken between lanuary and June with zenith angles 
ranging from 16 to 35 degrees. Analysis of these data with 
the standard MAGIC software (]Moralejo et al.|2009jl resulted 
in a 10 s.d. detection above 200 GeV. 

VERITAS — The Very Energeti c R adiation Imaging Tele- 
scope Array System (VERITAS j lu | consists of four 12m 
diameter IACTs and is located at the base camp of the 
Fred Lawrence Whipple Observatory in southern Arizona 
(31°40' N, 1 10°57' W; 1280 m above sea level). More details 
about VERITAS, the data calibration, and the analysis tech- 
niques can be found in Acciari et akj(j2008|. VERITAS ob- 
served M87 in 2010 for 48.2 h (after quality cuts) resulting 
in a detection of 26 s.d. The zenith angles of the observations 
ranged from 19 to 40°, with a few nights with zenith angles up 
to 60° during the flares (April 9-1 1th). The energy threshold 
of the analysis for the mean zenith angle of the observations of 
25° is 250 GeV. The VERITAS data covered the whole 2010 
flare period; a detailed study of the VHE spectral evolution 


109 http://www.mpi-hd.mpg.de/hfm/HESS/ 
1111 http://wwwmagic.mppmu.mpg.de/ 

1 1 1 http://veritas.sao.arizona.edu/ 


indicating a spectral change with flux will be published in a 
parallel paper ( Aliu et a k]201 1) >. 

2.2. High energy (HE) 

Fermi-LAT — The Large Area Telescope (LAT) onboard 
the Fermi satellite is a pair-conversion telescope that covers 
the energy range from 20 MeV to more than 300 GeV 
(HE) ( |Atwood et al.l|2009 >. The LAT instrument features a 
per-photon angular resolution of (/ 6S% = 0.8° at 1 GeV and a 
large field-of-view of 2.4 sr. The primary mode of operation 
is an all-sky survey mode, where the full sky is covered ap- 
proximately every three hours. The LAT data for this analysis 
consists of two years of nominal all-sky survey data between 
the energy range 100 MeV and 300 GeV, and spanning the 
mission elapsed time (MET) 239557417 to 302630530 
(August 4, 2008 through August 4, 2010). Event selections 
include “diffuse” class events recommended for point source 
analysis, a rocking angle cut of < 52°, and a zenith angle cut 
of < 100° in order to avoid contamination from the Earth’s 
limb. A 2 ks window beginning at MET 259459364 was also 
removed in order to avoid contamination from GRB 090323, 
which occurred nearby. All LAT analysis was performed 
using instrument response functions (IRFs) P6_V11_DIFFUSE 
and science tools v9r2®p®, along with the recommendec 112 
Galactic diffuse gll_iem_v®2_P6_VllT)IFFUSE. fit 
and corresponding isotropic spectral template 
isotropic_iem_v®2_P6_Vl 1_DIFFUSE . txt. 

An analysis of the LAT spectrum over the two-yea r period 
was performed using a binned likelihood analysis (Mattox 
et alT] 19961 selecting all events that fell within a 20° x 20° 
square region of interest (ROI) centered at the M 87 radio po- 
sition of the core. All point sources from an internal two-year 
preliminary catalog that fell within 15° of the source were 
included in the fit. All sources within the square ROI were 
modeled with a power law spectrum with their normalization 
and index as free parameters, while those that fell outside of 
the ROI were fixed to their catalog values. The M 87 spectrum 
was modeled as a power law with photon index and normal- 
ization parameters left free and using the radio position ( |Fey| 
et al.| |2004) as the source locati on. A point source was de- 
tected with a test statistic (TS; Mattox et al.||1996] > of 301, 
representing a detection of V301 ^ 17 s.d. From the result- 
ing fit, the photon index and flux (> 100 MeV) were found 
to be 2.16 + 0.07 and (2.66 + 0.36) x 10~ 8 phcm~ 2 s _1 , re- 
spectively. In order to test for curvature, the spectrum was 
also fit to a log parabola, where the TS of the overall fit was 
found to improve by only 0.89, which does not represent a 
statistically significant improvement over the single power 
law. The largest systematic errors can be attributed to un- 
certainties in the modeling of the diffuse background emis- 
sion. These were estimated by repeating the analysis with 
both binned and unbinned gtlike using a refined version 
of the current diffuse background model that is under devel- 
opment by the LAT collaboration. Systematic errors on the 
index and flux were thus found to be (+0.05/ - 0.01) and 
(+0.40/ - 0.13) x 10“ 8 phcm~ 2 s _1 , respectively. Compar- 
ing results from this analysis for the l ast 1 4 month with the 
initial 10-month spectrum reported in|Abdo et al.| (J2009]), no 
evidence of variability in the flux above 1 00 MeV was found. 
Comparing the flux above 1 GeV between these two epochs, 
however, a marginal indication of a rise in the flux of the latter 

1 12 http://fermi .gsfc.nasa.gov/ssc/data/access/lat/ 

BackgroundModels.html 
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epoch at a significance of 2s.d. is found. Systematic uncer- 
tainties on the LAT flux due to the instrument response are 
estimated at values of 10 %, 5 %, and 20 % above and below 
their nominal values at log(E/MeV) = 2, 2.75, and 4, respec- 
tively ([Atwood et al.|2009|. 

The LAT fight curve (> 1 GeV) was constmcted using the 
events that fell within a 10° circular ROI centered at the M 87 
radio position. Fluxes were derived for bins with a detec- 
tion significance exceeding 3 s.d., otherwise upper limits were 
calculated. To generate the fight curve, events were grouped 
into 56-day time bins, and a separate likelihood analysis us- 
ing gtlike was performed over each of the bins. All point 
sources from the two-year fit were included in the model over 
each interval. Sources that fell within the 10° ROI were fit 
with their normalization parameters free, while the photon in- 
dex of each source was fixed to the best-fit value obtained 
from the full two-year analysis. Both the index and normal- 
ization parameters of M 87 were left free, except in the case 
of upper limit calculations, in which case the photon index 
was fixed to the nominal two-year average of 2.16. In order 
to avoid modeling sources with a negative TS value, an initial 
fit over each interval was performed, and all sources found to 
have a TS< 1 were subsequently removed from the fit. Fol- 
lowing the method for variability detection outlined in |Abdo 
et al. (2010a i, the weighted average flux was first calculated 


with a resulting value of(1.62±0.18)xl 0~ 9 ph cm~ 2 s~ 1 . A x 1 
analysis was then performed by comparing the best-fit values 
of all points to the weighted average, and the resulting prob- 
ability P(x 2 > xlbs) was f° un d t° be 0.027, which represents 
a significance of 2.2 s.d. for the source to be variable and falls 
slightly below the threshold for variability defined in lAbdo 
[etak] ( |2010al >. 


2.3. X-ray 

Chandra — X-ray data have been taken with the Advanced 
CCD Imaging Spectrometer (ACIS) on board the Chandra 
satellite. For details of the Chandra data reduction proce- 


dures, see 

Harris et 

Harris et a 

. ( 

2009). 


used. This results in a frame time of 0.4s with 90 % efficiency. 


Although this setup was essentially free of pileup when Wil- 
lson & Yaiigl ( |2002 1 tested 


various options during 2000 July, 
with the advent of the ever increasing brightness of HST-1, 


pileup (Davis|2001| became a major problem so the measure 


of intensity was switched to a detector based unit: keV/s. This 
approach uses the event 1 file with no grade filtering (so as to 
recover all events affected by ’grade migration’) and energies 
from 0.2 to 17 keV are integrated so as to recover all the en- 
ergy of the piled events. Other uncertainties for piled events 
come from the on-board filtering, the ’eat-thy-neighbor’ ef- 
fect, and second order effects such as release of trapped charge 
(see Appendix A of Harris et al. 2009|). Although a small 


circular aperture was used for flux-map photometry in Har- 


ris et al.| (|2003|, the basic analysis for this pap er adopts the 


rectangular regions used in Harris et al. (2006) so as to en- 
compass more of the point spread function (PSF). All events 
within each rectangle are weighted by their energy and the 
sum of these energies, when divided by the exposure times, 
gives the final keV/s value used in the fight curve. Uncer- 
tainties are strictly statistical, based on the number of counts 
measured ( s[N /N) and typically range from 1% to 5%. 


2.4. Optical 

Hubble Space Telescope (HST) — Near Ultraviolet HST ob- 
servations were obtained with the Space Telescope Imaging 
Spectrograph (STIS) and the Advanced Camera for Surveys 
High Resolution Camera (ACS/HRC). STIS imaging was 
taken with the F25QTZ filt er that has its maximum through- 
put wavelength at 2364.8 A ( Kim Quijano et al. 2003 i. All ob- 
servations after August 2004 were taken with the ACS/HRC 
using two filters: F220W and F250W which have their max- 
imum throughput wavelength at 2255.5 A and 27 15.9 A, re- 


|2003j l. All science ready files 
I public 


spectively (Mack & Gilliland 
were retrieved from the STScl public archive and processed 
through the PYRAF task MULTIDRIZZEL (Fruchter et al. 
2009 ). Both HST detectors have a pixel scale of ~ 0.024" per 
pixel and a resolution that enables one to clearly separate the 
nucleus and the innermost components of the jet, i.e. HST-1. 
The details of these observations have been presented injPerl- 


|man et aL] ( |2003| ), |Madrid| ( |2009| ), and |Perlman et al.| ( |201 1 ) . 

Optical HST polarimetry observations were obtained with 
the ACS/HRC and the Wide-field Planetary Camera 2 
(WFPC2), using the F606W filter that has its maximum 
throughput at roughly ~6000 A and a nearly flat throughput 
curve from 4800 to 7000 A. The single WFPC2 observation 
(pixel scale ~ 0.1 arcsec/pixel) was obtained when the ACS 
was not operational, and used the POLQ polarizing filter, 
while the ACS/HRC observations used the POLVIS polariz- 
ing filter. For both of these polarizers, the F606W filter gives 
nearly optimal transmission of parallel polarized fight as well 
as rejection of perpendicularly polarized fight. Polarimetry 
observations were obtained at 18 epochs between December 
2002 and November 2007, with 14 of the 18 observations con- 
centrated between November 2004 and December 2006 on 
the same schedule as the ultraviolet photometry. As with the 
UV observations, all polarimetry was obtained from the HST 
archive and re-calibrated using updated flat field files and im- 
age distortion correction (IDC) tables ( Mobasher et al.|2002 


Pavlovsky et al.||2002). multidrizzle was used to combine 
and cosmic -ray reject the imag es, which were aligned using 
Tweakshifts (Fruchter et al.|2009), and CTE (Charge Transfer 
Efficiency) corrections were computed using data found in the 
ACS instrument handbook ( |Boffi et al.|20071 >. Once this was 
done, the images in each polarizer were combined according 
to recipes in the ACS and WFPC2 instrument handbooks, re- 
spectively (jBoffi et al.|2007||Biretta & McMaster|1997ji. Be- 
fore performing photometry and polarimetry, galaxy emission 
was subtracted from the images. This was done using a model 
computed in the Stokes I image using the IRAF tasks EL- 
LIPSE and BMODEL. Error bars for the polarimetry are typ- 
ically 2-3% for high signal-to-noise data (further details can 
be found in [Perlman et al.|2011j >. 

Liverpool Telescope (LT) — Hybrid R+V-band optical po- 
larimetry data (460-720 nm at FWHM) were taken with the 
2-meter Liverpool Telescope (LT j 113 [ located on La Palma, 
using the newly commissioned RING02 fast-readout imaging 
polarimeter (Steele et al .|2010) i. The polarimeter uses a rotat- 
ing polaroid with frequency of approximately 1 Hz, during the 
cycle of which 8 exposures of the source are obtained. These 
exposures are synchronized with the phase of the polaroid and 
following the analysis method of Clarke & Neumayer (|2002Jl 
allows determination of the degree and angle of polarization. 


http://telescope.livjm.ac.uk/ 
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The LT observations were taken during and shortly after 
the observed VHE high-state, MJD 55295-402, with three 
measurements taken contemporaneously to the time of the 
VHE flare, MJD 55295-97. Further data on M 87 have been 
taken since then and will be presented elsewhere. Total in- 
tegration times of typically 100 s were used in the observa- 
tions, corresponding to an achieved polarimetric accuracy of 
about 1 % for the brightness of the source. The field of view 
of the instrument is of 4 x 4arcmin with a pixel scale of 
0.45 arcsec/pixel, which at the distance of M 87 corresponds 
to a linear scale of ~40pc/pixel, equivalent to approximately 
half the distance between the core and the innermost jet com- 
ponent, the knot HST-1. Given the typical seeing during the 
observations of 1 .0 " to 1 .7 " (FWHM) the two components, 
core and HST-1, cannot be resolved individually. There- 
fore, an aperture radius for the integration of the signal of 
2.7 arcsec was used, so that both the nucleus and HST-1 were 
included. The outer jet is nevertheless well resolved and in- 
dependent light-curves were produced with the same aperture 
radius used for the core but now centered at knot-A, located 
12.3 arcsec downstream from the nucleus, revealing a steady 
polarization aligned with the jet direction. 

The greatest source of error in the determination of the po- 
larization levels is a systematic effect caused by contamina- 
tion by the bright host galaxy, which extends out to several 
kpc and dilutes the polarization signal from the nucleus. The 
strength of the contamination is estimated by adding a Sersic 
profile with n = 4 to a frame built by summing data from all 
Polaroid angles. This sh ows the well known flat tening in the 
central ~ 10 arcsec (e.g. Kormendy et a l.|2009 1 . Fitting the 
data from 9 to 2.7 arcsec reveals an excess within the central 
2.7 arcsec due to the core plus HST-1 of ~ 10 ± 5 %. The 
photometric flux of the background galaxy is of order 90 % of 
the total flux within the aperture used for extracting the polar- 
ization measurements. The measured fractional polarization 
was therefore multiplied by a factor of 10 and the flux mea- 
surement divided by a factor 10 in order to correct for this 
background light, and allow comparison with the high resolu- 
tion HST measurements of the core + HST-1. 

Light-curves for the polarization position angle (or electric 
vector position angle, EVPA, defined as 0.5 X arctanOi/g), 
where u and q are the linear Stokes parameters) were ob- 
tained for all the epochs of observation, and are presented in 
Fig. [4] along with the other polarization quantities. Due to 
host-galaxy contamination.the polarization position angle for 
the nucleus and HST-1 combined are measured to an accu- 
racy of 25°, considerably degraded in comparison to the 10° 
resolution achieved for the outer jet where the host galaxy is 
fainter. 


2.5. Radio 

Radio interferometers enable one to observe the jet of M 87 
with a large variety of angular resolutions and sensitivities, 
depending on the array size and observing frequency. In 
Tab. [T]the highest angular resolution achieved for each instru- 
ment contributing data to this paper is shown. Overall, they 
span an angular scale range from as large as ~ 0.1" (from 
the VLA at 22 GHz) down to a fraction of a milli-arcsecond 
(e.g. ~ 0.2 mas with the VLBA at 43 GHz, see Table 1). 
The longest baseline and highest frequency observations pro- 
vide the most valuable information about the compact, flat 
spectrum core; conversely, instruments with shorter baselines 
(such as the VLA) or lower observing frequency (like the 
EVN) are most valuable for the fainter and more extended 


TABLE 1 

Angular resolution for different radio observations. 


Instrument 

Wavelength 

Resolution 

VLBA 

43 GHz 

(0.21 x 0.43) mas 

VLBA (MOJAVE) 

15 GHz 

(0.6 X 1.3) mas 

VLBA 

2.3 GHz 

(7.5 X 3.9) mas 

VLBA 

1 .7 GHz 

(8.0 x3.4)mas 

EVN 

5 GHz 

(1.0 X 2.0) mas 

VLA 

15 GHz 

~ 0.13" X 0.12" 

VLA 

22 GHz 

~ 0.10" x 0.09" 


Note. — The beam sizes are given as full width at 
half maximum (FWHM). In case different resolutions 
have been used in the analysis the highest one is given 
in the table. 


HST-1 feature. 


VLBA 43 GHz — The 43 GHz VLBA (Napier et al. 1994) data 
from 2006 through 2008 were collected as part of an effort 
to study the d ynami cs of the inner jet near the launch region 
( |Walker et al. |2008) . The 2009 and 2010 observations were 
part of a project to find and study another VHE/radio flare 
like that seen in 2008 (Acciari et al. [2009| . The observations 
were made on the 10-antenna VLBA using a total bandwidth 
of 64 MHz. The data were reduced in AIPS following the 
usual procedures for VLBI data reduction including correc- 
tion for instrumental offsets using the autocorrelations, band- 
pass calibration based on strong calibrator observations, and 
correction for atmospheric opacity based on the system tem- 
perature data. The a priori amplitude calibration depended 
on the gains provided by VLBA operations, which are based 
on results from regular single-dish pointing observations of 
Jupiter, Saturn, and Venus averaged over many months. The 
images are based on data that are both amplitude and phase 
self-calibrated. The flux scale for each epoch was set by nor- 
malizing the self-calibration gain adjustments to the a priori 
gains for observations above 30° elevation on those antennas 
with good weather and instrumental conditions for that epoch. 
The flux densities typically accurate to within about 5%. 

Three VLBA flux densities are provided. The first is 
the peak brightness on the core in an image made with a 
0.21 x 0.43 mas beam. The second is the integrated flux den- 
sity within 1 .2 mas of the core and represents the total emis- 
sion from a region within a projected distance of 0.1 pc from 
the presumed position of the black hole at the radio core 
(0.1 pc = 340 R s for a 3.0 x 10 9 M G black hole). The third 
is the integrated flux density from the jet over the region 1.2 
to 5.3 mas from the core. 


VLBA 15 GHz ( MOJAVE ) — VLBA 15 GHz observations from 



VLBA 2.3 GHz — Observations at 2.3 GHz were made on 8 
and 18 April 2010 using 10 VLBA stations as a part of ex- 
periments BH163. Each session has a total on-source time of 
~15 minutes with the total bandwidth being 64 MHz. In order 
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to obtain better u, v coverage, the short scan blocks (~ 2 min- 
utes per block) were distributed uniformly at several hour an- 
gles. The initial data calibration was performed in AIPS based 
on the st anda rd VLBI data reduction procedures (the AIPS 
i. The amplitude calibration with the opacity cor- 


cookbool 
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rections was applied using the measured system noise tem- 
peratures and the elevation-gain curves of each antenna.The 
data were fringe-fitted and then averaged at short intervals 
(every 5 seconds and 1 MHz in time and frequency domains) 
before the imaging process in order to avoid the smearing ef- 
fects due to the time and bandwidth averaging at the HST-1 
region. The images were made in DIFMAP software ( |Shep- 


herd| 1997]> with the iterative phase/amplitude self-calibration 


processes. The off-source rms noises in the resultant images 
are 0.4-0. 5 mJy/beam. The peak flux densities are provided 
for the core region with the synthesized beam of 7.5 x 3.9 mas 
at -5°, and the integrated flux densities are provided for the 
HST-1 region. The errors in flux densities are assumed to be 
5% based on the typical VLBA calibration accuracy at this 
frequency. 


VLBA 1.7 GHz — 19 epochs of VLBA 1.7 GHz observa- 
tions from ~2005.0 - 2008.0 have been obtained in an ef- 
fort to monitor the evolution of the HST-1 knot following 


its brightening at X-ray, optical, and radio frequencies (Har- 
ris et al. 2009). The results of the first 9 of these observa- 


tions ( progra ms BH126, BH135) were presented in Cheung 
et al. 
the 


Ll <|2007l 
HST-1 co 


where apparent superluminal radio features in 
complex were discovered. For the additional 10 
observations (programs BC167, BH151), the identical cali- 
bration and imaging procedure was followed for this analysis. 
From maps restored with the same uniformly weighted beam 
(8.0 mas x 3.4 mas elongated north-south), we measured the 
peak core brightness (5% errors are assumed). Integrated flux 
densities (10% errors are assumed) for the HST-1 were mea- 
sured from naturally weighted images (11.5 mas x 5.5 mas) 
using a 80 mas x 50 mas box that covers the entire radio com- 
plex resolved in the VLBA images (cf., Cheung et al. 2007)>. 


EVN — M87 was observed with the European VLBI Net- 
work (EVN) at 5 GHz seven times between November 2009 
and lune 2010 as part of a project aimed at correlating the 
radio and high energy behavior of the source. Given the inter- 
esting episodes of activity of the source during the campaign, 
a few observations were scheduled as ToO and do not have 
the same array configuration as the other ones. In general, 
6 to 11 telescopes participated in the observations, with base- 
lines ranging between less than 100 km (as provided by MER- 
LIN stations) and 7 000-9 000 km (as provided by the Arecibo 
and Shanghai stations, respectively). The observations were 
carried out at 5.013 GHz, divided in 8 sub-bands separated 
by 16 MHz each, for an aggregate bit rate of 1 Gbps. The 
data were correlated in real time at JIVE using the so-called 
e-VLBI technique, which provides a fast turn-around of the 
results (as was the case near the Feb. 2010 event; see |Giro-] 
letti et al .|2010| >. Automated data flagging and initial ampli- 
tude and phase calibration were also carried out at JIVE us- 
ing dedicated pipeline scripts. Data were finally averaged in 
frequency within each IF, but individual IFs were kept sep- 
arate to avoid bandwidth smearing. Similarly, the data were 
time-averaged only to 8 s, in order to avoid time smearing. 
Final images were produced in DIFMAP after several cycles 
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TABLE 2 

Parameters of the fit to the 2010 VHE flare 
data 


Parameter 

Value 

Unit 

Fit range 

55290 - 55299 

MJD 

-.rise 

1.69 + 0.30 

d 

decay 

T d 

0.611+0.080 

d 

to 

55295.954 ± 0.094 

MJD 

®0 

(2.01 ±0.15) x 10 -11 

-2 -1 

cm s 


Note. — The fit results in a ^/d.o.f. = 
10.02/11 with chance probability P = 0.53. The 
parameters are defined in the text. 


of phase and amplitude self-calibration. Various weighting 
schemes were applied to the data to improve resolution in the 
core region and enhance the fainter emission in the HST-1 
region. For the core, uniform weights were used obtaining a 
typical HPBW of 1. 0x2.0 mas (in PA -20°) and providing the 
peak brightness. For HST-1, natural weights were used result- 
ing in a 7.5 x 8.5 mas HPBW; a resolved structure is clearly 
detected and integrated flux densities for the full region are 
given. For additional details, please see Giroletti et al. (201 1[). 


VIA — Data from the VLA archive for 10 epochs have 
been used, selecting observations performed in A-array at 15 
and/or 22 GHz. The typical angular resolution of the VLA 
is ~ 0.13" x 0.12" and - 0.10" x 0.09" at the two frequen- 
cies, respectively, which permits one to resolve the core and 
the HST-1 feature. Data reduction was carried out in AIPS 
in the standard manner: the flux density scale was tied to the 
main gain calibrator 3C 286 (which was used in all the ob- 
servations) using SETJY with the available model. Phases 
were calibrated using the compact source 1224+035, obtain- 
ing good solutions. A few final iterations of the phase and 
amplitude self-calibration have been carried out to improve 
the image quality. Integrated flux densities for the core and 
HST-1 have been derived using JMFIT. 

Additional 15 GHz VLA data for HST-1 from lHarris et al.l 
(2009) are also shown in the light curve. 


3. THE 2010 VHE CAMPAIGN 

The 2010 flare — During the 2010 VHE monitoring campaign 
two episodes of increased VHE activity have been reported 
( Mariotti & MAGIC Collaboration |20 1 0b[ Ong & Mariotti| 
|2010| i: The first episode took place in Feb. 2010 where a sin- 
gle night of increased activity was detected by MAGIC (Fig. [2] 
bottom panel around -35 d; detection significance > 5s.d.). 
Follow-up observations by H.E.S.S. and VERITAS did not 
reveal further activity at VHE. The second episode took place 
in Apr. 2010 and showed a pronounced VHE flare detected by 
several instruments triggering further MWL observations. In 
the following, the discussion will concentrate on this second 
flaring episode. 

The VHE activity of this second flaring episode is concen- 
trated in a single observation period between MJD 55290 and 
MJD 55305 (-15 days; see Fig. [2] bottom panel & Fig. [3]). 
This time -period is exceptionally well covered with 21 point- 
ings by different VHE instruments, resulting in an observation 
almost every night. It should be noted that during nights with 
(quasi) simultaneous observations by different instruments, 
the measured fluxes are found to be in very good agreement. 

The detected flare displays a smooth rise and decay in flux 
with a peak around MJD 55296 (April 9-10, 2010). A peak 
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Fig. 2. — VHE light curve of M 87 of the flaring episodes in 2005 (top), 2008 (middle), and 2010 (bottom). Integral fluxes are given above an energy of 350 GeV. 
The lengths of the gray bars correspond to the length of the gray shaded areas in Fig.|T] A time of 0 days corresponds to MJD 53460, MJD 54500, and MJD 55270 
for 2005, 2008, and 2010, respectively. Flux error bars denote the 1 s.d. statistical error. Horizontal error bars denote the time span the flux has been averaged 
over. Note that in the case of time spans longer than one night the coverage is not continuous. 


flux of ~ 2.5 x ICE 11 cm _2 s _I is reached, which is about a 
factor 10 above the quiescent flux level of the source. The 
data on the rising part of the flare indicates a steady rise. On 
the decaying side the situation is more complex: two nights 
after the detection of the maximum flux (MJD 55298) all three 
instruments measured a low flux compatible with zero, while 
in the following 3 nights (MJD 55299/55300/55301) a higher 
flux is detected by VERITAS. 

To derive the timescales of the flare the VHE light curve 
from MJD 55290 to MJD 55299 is fitted with a two-sided 
exponential function 


<p = cp 0 x e ,oI/At with 


( At = Ar rise for t < t 0 
\At = Ar decay for t > t 0 


(1) 


resulting in an excellent description of the data with a 
^f 2 /d.o.f. = 10.02/11 and a corresponding chance probabil- 
ity of P — 0.53. The resulting fit parameters are summarized 
in Tab. [2] Very similar results are obtained when (i) fitting 
a longer time span (e.g. MJD 55275 - 55305) or when (ii) 
adding a constant offset as a free parameter to account for 


a possible baseline flux, though the chance probabilities of 
these fits are slightly worse then for the original fit described 
above. 

Flux doubling times of Td = ln( 2) x At of r™ e = (1.69 + 

0.30) days for the rising and rj ecay - (0.61 1 ± 0.080) days for 
the decaying part of the flare are derived, which significantly 
differ by a factor 2.77 ± 0.62. The variability timescales de- 
rived for the 2010 VHE flare are the most precisely measured 
VHE variability timescales determined for M 87 to date. Pre- 
viously detected VHE flares only allowed for rough estimates 
due to the variability pattern, sampling, and statistics. 

Comparison with the 2005 and 2008 flares — The VHE light 
curve from 2005, 2008, and 2010 around the flaring episodes 
is displayed in Fig. [2] In all three flares, similar flux levels are 
reached. The apparent timescales of the order of days are also 
very similar. The time period over which activity is detected, 
is comparable in 2008 and 2010 but is slightly longer in 2005, 
though the gaps in the sampling make any definitive statement 
difficult. During the 2010 activity period a pronounced flare is 
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Fig. 3. — VHE light curve of M 87 zoomed on the 2010flare. Also shown are 
the results of the fit of an exponential function to the data. Error bars denote 
the 1 s.d. statistical error. The results of the fit to the data are summarized in 
Tab.[2] Further details can be found in the text. 


detected which is well described by an exponential behavior 
(see previous paragraph). While the 2005 flare is compatible 
with such behavior (the sampling is considerably worse than 
in 2008 and 2010), the 2008 activity state seems more erratic 
with several maxima over a similar time period as in 2010. 
Taking the best fit function to the 2010 flare as a template, 
large parts of the 2008 flaring activity are not compatible with 
the 2010 behavior (see the Appendix for details). 


VHE flare duty cycle — Given the typical length of a VHE 
flare of order a day, the nightly VHE y-ray flux measurements 
presented in this paper can be used to estimate the duty cycle 
of M 87 for VHE flares^ For exa mple , following the approach 
presented in Jorstad et al. (2001 one derives a duty cycle 
of ~28%. However, due to the observing strategy followed for 
M 87 (observations have been intensified after the detection 
of a high state), the data-set is biased. The derived duty cycle 
is, therefore, likely overestimated and should be considered 
an upper limit. To enable a more unbiased view and to give 
an estimate of the possible uncertainties the duty cycle for 
a range of threshold fluxes is presented, calculated from the 
number of data points above the threshold flux divided by the 
total number. For a threshold flux of ^thresh = (0.5/0. 8/1) x 
10~ H cm~ 2 s -1 the respective duty cycles are ~14% / 7% / 
4%. 


4. OPTICAL POLARIMETRY 

Radio and optical emission from regions within the M 87 
jet show a high degree of polarization (Perlman et al. ! 1 999| >, 
which is indicative of non-thermal emission processes in 
highly ordered magnetic fields (synchrotron radiation). In 
general, polarization and especially changes in the polariza- 
tion can be utilized to investigate the magnetic field struc- 
ture of the emission sites and put additional constraints on 
the emission models. Recently, changes in the polarization 
(amplitude and angle) have been used to investigate the origin 


of the y-ray emission in several blazars (Marscher et al. 

2008. 

[2010; Abdo et al . |20 1 0b ; Barres de Almeida et al/pOlO 

). The 


116 Calculating the number of observations with a flux a factor 1.5 above 
the error weighted mean (for this data set ~ 1.8 X 10 1 - cm " 2 s 1 ). 


observed pattern - a change in the HE/VHE flux coinciding 
with a rotation in the electric vector position angle (EVPA) 
and/or a change in the relative polarization - could result 
from: a non-axisymmetric magnetic field through which the 
emission region propagates, a swing of the jet across the line 
of sight, or a curved trajectory of the emission region e.g. fol- 
lowing helical magnetic field lines or a bent jet (lAbdo et al. 
12010bD . 


HST 2001-2008 — Regular optical polarimetry observations 
of the M 87 j et have been performed w ith the HST from 2004 
to 2006 (Sec.|2.4| Perlman et al.|201 1|). The measured degree 
of polarization and the polarization angle for the core region 
and HST-1 are shown in Fig. [4} left panel. While both compo- 
nents show variable polarization their overall behavior is very 
different: during the flaring period (2004-2007) HST-1 shows 
a clear correlation of the fractional polarization with the total 
flux, with the fractional polarization ranging from 20 to 40 %. 
The EVPA remains almost constant ~ 7 - 8° off the direction 
of the jet (-69.5°; gray line in Fig. [4] lower panel). The nu- 
cleus, on the other hand, displays a highly variable EVPA with 
a lower overall polarization of 2 to 10 %. The EVPA changes 
are not correlated with the total flux or the fractional polariza- 
tion and appear, within the given sampling, erratic. The HST 
polarimetry data covers the 2005 VHE flare but, apart from 
the previously mentioned correlation of the HST-1 photomet- 
ric flux and fractional polarization, no correlations connected 
to the 2005 VHE flaring episode are apparent from the data. 


LT2010 — Triggered by the detection of the 2010 VHE flare 
the LT started to take regular optical polarimetry observations 
of M87 in April 2010. The resolution of the instrument is 
not sufficient to separate the core and HST-1 and, therefore, 
one value for the combined core + HST-1 region is given. 
The LT detects clear variability of the degree of polarization 
of the core + HST-1 complex and marginal evidence (~ 2cr) 
for variability of the EVPA (Fig. |4l right panel). The polar- 
ization degree changes between ~ 2 and ~ 12% but without 
any apparent correlation with the VHE y-ray flux. The mea- 
sured EVPA indicates changes at the time of the VHE flaring 

episode from 130° to ~ -60° (~ 70° in 3 days) and later 

settles to a stable EVPA of ~ -20° about 50 days after the 
flare, though the only marginal indication for EVPA variabil- 
ity and the lack of continuous monitoring before and directly 
after the VHE flare makes it difficult to interpret these changes 
in terms of a continuous rotation of the EVPA. The observed 
polarization variability pattern is compatible with the pattern 
previously observed by HST for the core region and different 
to what has been observed from HST-1. This indicates that 
the 2010 LT data is possibly dominated by emission from the 
core. 


5. VHE FLARES & MWL CORRELATIONS 

Over the last 10 years M 87 has been extensively monitored 
all across the electro-magnetic spectrum from radio to VHE 
(Fig. [TV This large data-set can be used to investigate MWL 
correlations and thereby probe the origin of the VHE y-ray 
emission. In principle, many different physics processes could 
contribute to the production of VHE y-ray emission in M 87, 
e.g. annihilation of massive dark matter particles (jBaltz et al.| 
2000 ), cosmic ray interactions in the extended radio lobes and 
the surrounding cluster (jPfrommer & EnBlin ||2003 jl , particle 
acceleration in the relativistic jets, etc. The detected short- 
term variability with timescales of the order of days and the 
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Fig. 4. — Optical polarization for different components of the M87 jet as measured by HST and LT. The LT p olari zation measurement of the core + HST-1 
complex is scaled up by a factor of 10 to compensate for the host galaxy emission in the measurements (see Sect. |2.4| for details). The horizontal gray solid line 
in the lower panel marks the jet direction (—69.5°) while the horizontal gray dashed lines are located at -69.5° ± 9U U (180° distance). 


limits on the location place strong constraints on possible sce- 
narios (Aharonian et al. 2006a} , leaving the inner jet and the 
close vicinity of the SMBH as the most probable emission 
sites. In the following, the MWL behavior of the two most 
prominent features in the innermost structure of M 87, namely 
the HST-1 knot and the core, are discussed in the light of the 
VHE flaring activity. 


HST-1 — Between 2001 and 2008 the first bright feature in 
the jet resolved by the HST in the optical, HST-1, underwent 
a spectacular flare detected in radio , optical, and X-rays (Perl- 
man et al.|2003||Harris et al.|2003}. The flare displayed a rel- 
atively smooth rise over several years with a flux increase by 
more than a factor 50 in X-rays and optical. The flux peaked 
in the beginning of 2005 (Hams et al.| 2006| 2009} around 
the same time when the enhanced activity level a nd the first 
short term v ariability had been detected at VHE (Aharonian 
et al. |2006a] ). HST-1 has been discussed as a p ossible site 
for the VHE y-ray emission (e.g. Stawarz et al.| 2006 |Che 


ung et aL|2007| Harris et al.||2009 1 . While the size of 
1 as a whole is too large to account for the short-term vari- 
ability detected at VHE (following causality arguments), high 
resolution VLB A radio observations resolve HST-1 into sev- 
eral, partially unresolved sub-structures (Cheung et al.|2007}. 
These sub-structures also display apparen t sup erluminal mo- 
tion up to 4.3 c ± 0.7 c (see also Giovannini et al~| [20 1 1 } . In 
combination with the detected synchrotron X-ray emission 
and strong polarization of the radio-to-optical continuum, this 
indicates that efficient in-situ acceleration of the radiating par- 
ticles is taking place in compact sub-volumes of the HST-1 re- 
gion, characterized by well-organized magnetic field and rel- 
ativistic bulk velocities. On the other hand, during the 2008 
and 2010 VHE flares HST-1 was in a low flux state without 
pronounced activity at radio or X-ray wavelengths, thus disfa- 


voring it as the origin of the VHE y-ray emission during these 
episodes. 


Core — The direct vicinity of the SMBH and the jet base 


y-ray emission (e.g. 

Reimer et al. 2004; Ghisellini et al. 

2005 

Tavecchio & Ghise 

lini 2008[]Neronov & Aharonian 

2007 

Rieger & Aharonian 

2008, Lenain et al. 

2008, Barkov et al. 

2010; Giannios et al. 2010 Levinson & 

Riege 

r 201 1 ). 

M87 
et al. 

is only a weak IR source ( vL v ~ 10 41 erg/s; 

Perlman 


cape even from the close vicinity of the SMBH without suf- 
fering strong absorption due to yy-interactions (Neronov & 


Aharonian|2007||Brodatzki et al.|201 lb, although the debated 
origin of the observed IR photons and the poorly k 


origin of the observed IR photons and the poorly known struc- 
ture of the accretion disk in the M 87 core, which both play a 
crucial role in this context, should be kept in mind (see the dis- 


cussion in Cheung et al. 

2007 and Li et al. 2009). The M87 

jet base has been imagec 

with VLBI with sub-mas resolution 


resolved, edge brightened structure to within 0.5 mas of the 
core (0.04 pc) and indications for a weak counter-jet feature, 
suggesting that the SMBH lies within the central beam of the 
VLBI observations. At even shorter distances to the core the 
jet has a wider opening angle, which is interpreted as the jet 
collimation zone ( Junor et al. 1999} . 

In 2008, densely sampled 43 GHz radio observation of the 
innermost jet regions revealed a flare of the radio core (flux 
increase of ~ 30%;|Acciari et al. |2009} . At the same time, a 
flare at VHE and a subsequent enhanced X-ray flux from the 
core region followed by a sharp decrease were detected. The 
observed MWL variability pattern supported the interpreta- 
tion that the VHE y-ray emission likely originates from the 
close vicinity of the SMBH near the jet base (Acciari et al. 
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[2009). The observed MWL behavior (VHE and radio flux, 
radio map) is well described by a simple, phenomenological 
model where the VHE flares are associated with the injection 
of plasma at the jet base (Acciari et al7| [2009[ supporting on- 
line material). As the injected plasma blobs travel down the 
jet, they expand and become transparent for radio emission 
leading to the observed radio feature. 

In contrast, radio observations taken in 2010 contempora- 
neous with the VHE flare show no enhanced radio flux from 
the core region (Fig. [TJ: VLBA 43 GHz ToO observations 
triggered by the detection of the VHE flare indicate a sta- 
ble flux state of the core (1.2mas) and the inner jet (1.2 to 
5.3 mas) at the previously detected flux levels. In addition, 
EVN 5 GHz, VLBA 2.3 GHz, and MOJAVE 15 GHz measure- 
ments also show no indication of an enhanced radio flux state 
from the core in 2010. A direct comparison of these data with 
the 43 GHz data is difficult given the differenc e in resolution 

In the op- 
and HST-1 


and the missing overlap during the 2008 flare 
tical, only a combined measurement for the core 
with limited sensitivity by the LT is available, which does not 
indicate any strong activity in the two components in 2010. 

On the other hand, Chandra X-ray observations of the core 
show an enhanced flux ~ 3 days after the peak of the VHE 
y-ray emission in 2010 (see Fig. [I). The flux is enhanced by a 
factor ~ 2 for a single measurement and than drops back to a 
lower state less than two days later. The observed variability 
timescale is significantly shorter (by a factor ~ 10) then the 
shortest X-ray variability measured previously from the M 87 
core (< 20 days; |Harris et al.] 2009[ ). Further details on the 
Chandra X-ray data from 20 10 will be reported in a compan- 
ion paper (Harris et al.[201 1 1 . 

It should be noted that the X-ray fluxes measured from the 
core during the time of the VHE flares in 2008 and 2010 are 
the two highest measurements since the start of the Chandra 
observations in 2002. This coincidence can be interpreted 
as indication that the VHE y-ray emission in 2008 and 2010 
originates from the X-ray core region. During the 2005 VHE 
flaring episode no enhanced X-ray emission from the core 
was detected. At that time, HST-1 was more than a factor 30 
brighter than the core region in X-rays leading to uncertain- 
ties in the flux estimation of the core (e.g. ’eat-thy-neighbor’ 
effect; see Appendix A of Harris et al.|2009| which could sup- 
press the detection of a core flare. 

6. DISCUSSION 

General considerations — The VHE and MWL data in this 
paper can be interpreted in two fundamentally different ways: 
(1) each of the VHE flares detected originates from a different 
emission region and/or process, or (2) they have a common 
origin. Support for the former interpretation comes from the 
difference in the VHE light curves (variability pattern and the 
overall duration of each episode), as well as from the differ- 
ence in the apparent MWL correlations in the different years. 
In this interpretation, all previously derived models remain 
possible and an additional explanation for the 2010 flare has 
to be found. Observational support for the latter interpreta- 
tion, involving a common origin of the VHE flares, follows 
from the similarities in the detected peak fluxes, flux doubling 
timescales, and also in the spectral shapes of the different 
VHE flares (for details on the VHE spectrum see |Aharonian 
|et al.|2006a[| Albert et al.|2008a[|Aliu et al.|2011) . In this in- 

117 Noteworthy, one of the highest flux states in the MOJAVE 15 GHz data 
is measured in 2008 shortly after the VLBA 43 GHz observations ended. 


terpretation the only remaining MWL signature for the VHE 
flares is an X-ray flare of the core region (see previous section) 
locating the VHE y-ray emission site in the central r esol ution 
element of the Chandra observations (0.6 " ~ 50 pc) 
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In addition, it is not known whether the quiescent and flar- 
ing VHE y-ray emission are produced in the same region and 
by the same process. While, in principle, they can be pro- 
duced in two distinct locations by two different processes, the 
similarities of the VHE y-ray spectrum between the two states 
might suggest a common origin, or at least a very similar 
physical process involved. A common origin of the quiescent 
and flaring VHE fluxes is therefore anticipated below, even 
though different emission regions dominating the two states 
cannot be excluded. Interestingly, the HE spectrum measured 
by Fermi-LXI joins smoothly with the VHE spectra derived 
for the quiescence state (jAbdo et al.|2009) possibly indicating 
a common origin and, therefore, HE-VHE flux correlations 
could be expected. Unfortunately, due to the low flux of M 87 
in the HE band and the limited sensitivity of Fermi-LAl, no 
conclusive statement on such correlation can be derived from 
the data presented in this paper. 

In the following, the discussion will mainly focus on the 
case where a common origin of all VHE flares in the X-ray 
core (0.6 " ~ 50 pc) associated with an X-ray flare is assumed. 

Characteristics of the VHE flares — The VHE flares are char- 
acterized by a variability timescale f var - 1 day, and a broad 
band power law spectrum extending from ~0. 1 Te V up to 
at least lOTeV with photon index T - 2.3, weakly vari- 


able during flares 

Aharonian et al. 2006a, Albert et al. 

2008a| Aliu et al. 2 

Oil 

). The observed VHE flux during 
ted similar maximum flux levels of 

the haring episodes reac 


®>0.35TeV - (1 - 3) x 10 11 ph cm 


With the adopted 


distance of d — 16.7 Mpc and using the average photon in- 
dex T = 2.3, the corresponding isotropic VHE luminosity is 
Lvhe = 47rof 2 a> >£o £o(r-l)/(r-2) - (0.8-2.4)xl0 42 ergs-'. 
If extrapolated down to the HE range (0. 1 Ge V) with the 
same photon index, the total y-ray luminosity of the flaring 
events would be even higher, namely Lhe-vhe ~ 10 43 erg s 1 . 
This is a non-negligible amount given that the bolometric 
accretion luminosity of the M 87 nucleus is relatively low. 


10 42 ergs 1 


Di Matteo et al. 


(Reynolds et al. 1996 
1 2 00 3 [ |Kharb & Shastri||2004[ ), and the total kinetic lu mirioi 
ity of the j et is also quit e modest, L\ ~ 1 0 44 erg s 1 (Bick- 
[nell & Begelmann| 1 996; Owen et al. 2000) . The efficiency 
of the VHE y-ray production in M87 is therefore an issue, 
even taking into account order-of-magnitude dimmer quies- 
cence fluxes. More specifically, the observed VHE luminos- 
ity during the flaring events (which is approximately equal to 
the average/quiescence bolometric y-ray luminosity in both 
HE and VHE bands; see |Abdo et al.| 2009] > is of the order 
of the accretion luminosity in the M 87 system, constituting 
at the same time about 1% of the jet total kinetic luminosity, 
Tvhe ~ T acc — 0.01 x Lj. 

Interestingly, assuming the observed emission is moder- 
ately beamed with a Doppler factor 6 of the order of a few 
and that, when viewed at smaller inclinations, the beaming of 
the nuclear jet in the M 87 system would be the same as the 
one deduced for blazar sources, namely 6* - 10-30, the 
isotropic VHE flaring luminosity of M 87 observed at smaller 
viewing angles would read L* HE — (d*/d) 4 Lv he — ( 10 44 - 

118 It should also be noted that direct observational limits on the extent of 
the VHE emitting regio n cons train its projected size to be <, 14kpc centered 
on the radio core l|Aharonian et al.|2006aj. 
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10 4( ’)ergs This is in the range of VHE flaring luminosi- 
ties of blazars of the BL Lac type, for which low -power radio 
galaxies like M 87 are believed to constitute a parent popula- 
tion ( jUrry & Padovani||l995) l. Moreover, in such a scenario 
the observed variability timescales are scaled down by the ra- 
tio of the Doppler factors At* = (6/6*) At ~ (0.5 - 0.2) A t 
becoming of the order of a few hours, which is less than the 
timescale derived from the size of the Schwarzschild radii 
via causality arguments (-day). Variability timescales shorter 
then the causality timescales are also observed in flaring VHE 
blazars, where, in some cases, variability timescales down to 
a few minutes imply even higher Doppler factors of 0(100) 


T ; "Lvhe, 


(e.g. Albert et al. 2007] |Wagner|2008| |H,E S.S Collaboration 
et al.|2010[ ). This agreement can be considered as support for 
a blazar-like origin of the y-ray emission in M 87. The VHE 
flux changes by a factor of 10 and a VHE flare duty cycle 
of 0(10%), as estimated in this pape r, would also be consis- 
tent with such a hypothesis (Wagner|2008|l. Note in this con- 
text that with relativistic beaming involved the power emitted 
during the VHE flares would be reduced with respect to the 
isotropic luminosity estimated above as L em , vhe ~ r ~ 2 
where V j is the bulk Lorentz factor of the emitting region. 

Particularly interesting and constraining are the distinct 
timing properties of the one-day-long VHE flares of M 87, 
recognized clearly and quantified for the first time in this pa- 
per. These include repetitive outbursts (2005), erratic flux 
changes (2008), exponential flux increases/decays for well- 
defined isolated events, and significantly different rise and de- 
cay times (2010). If indeed the three flares share a common 
origin, every model dealing with the VHE y-ray emission of 
M87 has to be able to accommodate such a diversity in the 
observed variability patterns. 

Lastly, MWL correlations (or the apparent lack of such) 
have to be addressed as well, meaning in particular the emerg- 
ing connection between the VHE and X-ray bands with no ac- 
companying flux variations at radio or optical frequencies. It 
should be noted, however, that unlike the order-of-magnitude 
flux increases observed in the VHE regime, the observed vari- 
ability in the X-ray regime is of a much lower amplitude 
(flux changes by a factor of ~ 2 only), corresponding to a 
rather moderate X-ray core peak luminosity of the order of 
~ 10 41 erg s 1 . This value has to be taken with some caution 
though since, up to now, there are no truly simultaneous X-ray 
core observations during a VHE flare (i.e. during the night of 
the peak of the VHE y-ray emission). 

In the following paragraphs, existing models for VHE y- 
ray emission from the M87 core are briefly reviewed and 
discussed in the light of the new observational result pre- 
sented in this paper. The majority of models discussed have 
been published before the detection of M 87 by Fermi-LAT 
( Abdo et al.||2009 1 . Models published before 2006 are also 
not constrai ned by the VHE short-te rm variability discovered 
by H.E.S.S. (Aharonian et al .[2006a| . 


Observations versus modeling — Let us start with the ’mis- 
aligned blazar’-type, or rather ’misaligned BL Lac’ -type sce- 
narios (noting at the same time that the standard leptonic ’ho- 
mogeneous one-zone SSC’ approach has been often consid- 
ered to fail in explaining the observed VHE properties of the 
M 87 nucleus; but see also the discussion in |Abdo et al.|2009jl. 
One of the first models of this kind was discussed by Reimer 
et al. (2004}l, and involved a mixture of hadronic and leptonic 
emission processes operating within its relativistic outflow at 
distances from the jet base small enough to assure R < c f val 6 


for the emission region size R ~ 10 16 cm, variability timescale 
t v ar ~ 1 day, and the expected moderate beaming, 6 ~a few. 
Designed to account for the non-simultaneous data known at 
that time, the HE/VHE continuum was dominated by the syn- 
chrotron radiation of protons and muons in strong magnetic 
fields. While proton synchrotron radiation, with a maximum 
intrinsic cut-off at ~ 0.3 TeV, appears unable to account for 
the measured VHE spectra detected up to ~ 10 TeV, syn- 
chrotron radiation from muons (produced in charged pion de- 
cays) and the pion cascade components could potentially ex- 
tend the VHE spectrum to higher energies. This would then 
imply pion-production losses to be at least of the same order 
of magnitude as proton synchrotron losses. The low energy, 
millimeter-to-UV part of the spectrum is ascribed to the syn- 
chrotron emission of primary electrons with the same injec- 
tion spectral index as the primary protons, and the high, non- 
simultaneously measured X-ray flux considered to be domi- 
nated by an additional component. A direct VHE/X-ray cor- 
relation with a high X-ray flux level is thus not expected in 
such a sc enario. 

Georganopoulos et al. (2005j> discussed a leptonic version 


of the blazar-type modeling of the M 87 nucleus, relaxing the 
assumption regarding the homogeneity of the emission re- 
gion. More precisely, [Georganopoulos et al. considered the 
case of a relativistic jet decelerating substantially on sub-pc 
distances from the core, and showed that the velocity dif- 
ference between a faster and a slower portion of the outflow 
could lead to the enhanced inverse-Compton emission of the 
former one in the TeV range (due to a relativistic boost of the 
energy density of soft photons produced in one portion of the 
flow in the rest frame of the other), at a level allowing a fit 
to the VHE data. In the model, the core emission observed at 
longer wavelengths (including the GeV range) was predomi- 
nantly due to the slower and outer parts of the decelerating jet. 
A related scenario was analyzed by Tavecchio & Ghisellini 
(2008). who assumed instead radial velocity stratification, in- 
volving a relativistic jet spine surrounded by a slower sheath. 


In the particular model fits presented by Tavecchio & Ghis 
|ellini| the observed radio-to-HE emission of the unresolved 
M 87 core was dominated by the radiative output of the jet 
spine, while the observed VHE y-ray emission was produced 
mainly in the jet boundary layer. Both models are charac- 
terized by compact sizes of the VHE y-ray emission region 
(R < 10 17 cm), moderate beaming, and no obvious — if any 
— correlation between the VHE and lower-frequency bands. 
Importantly, the VHE spectra calculated in the framework of 
both models were rather soft, due to the Klein-Nishina and y- 
ray op acity effects involved (see the discussion in [Tav ecchio 


& Ghisellini|2008 1 . Multi-zone approaches, meaning a dou- 


bled number of the model free parameters, could possibly al 
low for some modifications to the presented fits enabling one 
to accommodate some additional MWL constraints (for ex- 
ample the VHE/X-ray correlation), but the observed flat flar- 
ing spectra of M 87 in the VHE range will remain problematic 
in both cases. 

Blazar type models assuming a different s tructure of the 
emission region are discussed in Lenain et al.|(|2008|) and|Gi- 
annios et al.| ( |2010| >. In those, the bulk of the observed emis- 
sion of the unresolved M 87 jet was proposed to be produced 
in extremely compact sub-volumes of the main outflow. In the 
scenario analyzed by Lenain et al. such compact sub-volumes 
were assumed to be multiple blobs ejected from the core with 
modest or highly relativistic bulk velocities into a cone with 
a large opening angle. In the framework of the model con- 
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sidered by Giannios et al. on the other hand, one is dealing 
with fast mini-jets moving in random directions within a rel- 
ativistic ‘large-scale’ outflow. Complex setups of the models 
resulted in the fact that different variability patterns and vari- 
ous variability timescales could, in principle, be expected, de- 
pending on a particular choice of the model free parameters 
(for which broad ranges of values could be considered). In the 
particular fits presented by Lenain et al7| (|2008|>, the authors 
considered for example the radio-to-UV continuum of the 
M 87 core to be produced by the extended though still unre- 
solved portion of the jet, and argued that the synchrotron and 
inverse-Compton emission of tiny blobs moving within such 
a jet may account for the observed nuclear X-ray and y-ray 
fluxes, respectively, including not only the quiescent but also 
flaring states. The anticipated VHE/X-ray correlation with no 
accompanying flux variations at longer wavelengths is thus an 
interesting feature of the model. 

Yet the crucial assumption regarding linear sizes of the 
emitting blobs significantly smaller than the Schwarzschild 
radius of the M 87 black hole, R <s Rs, might be regarded 
as questionable (see the discussion in Begelman et al. 2008| >. 
A possible physical justification for the formation oF"such 
ultra-compact sites of the enhanced energy dissipation close 
to SMBHs in AGN jets, and M 87 in particular, was given by 
Giannios et al. (j2010|, who speculated that this is the mag- 
netic reconnection process which may trigger compact beams 
of plasma (‘mini-jets’), which then propagate with relativis- 
tic bulk velocities within a strongly magnetized, extended 
and similarly relativisti c outflow. In their ap plication to the 
specific case of M 87, |Giannios et akj (j2010|! demonstrated 
that the synchrotron and inverse-Compton emission of the 
reconnection-driven mini-jets could account for the observed 
X-ray and VHE nuclear fluxes, somewhat in analogy to the 
results obtained by Lenain et al7| ( |2008) >, again under particu- 
lar model assumptions regarding the model free parameters. 
In the end, not only the fast VHE variability and the flat VHE 
spectra, but also the VHE/X-ray correlation could be success- 
fully accommodated in the model. 

Possibly challenging for scenarios involving ultra-compact 
emission regions is also the required efficiency of the VHE y- 
ray production. As discussed above, the VHE flares observed 
in the M 87 system are characterized by an isotropic luminos- 
ity constituting at least 1% of the total kinetic power of the 
M 87 jet. Whether such a power can indeed be efficiently dis- 
sipated within tiny sub-volumes of the jet remains an open 
question, even taking into account the expected beaming cor- 
rections. 

Magnetospheric particle acceleration and emission models 
have also been invoked to explain the observed VHE emis- 
sion in M87 (see e.g. Rieger [20 1 1 1 . for a review). Assumed 
to be operating close to the event horizon of the central su- 
permassive black hole, the anticipated variability timescale 
(light-crossing argument) could be as small t - 2G;V/bh/c 3 - 
0.35 day, and thereby satisfy the observed variability con- 
straints. The observed radio-VHE correlation during the 2008 
VHE flare has been interpreted to provide additional support 
for such an approach (Acciari et a l.|2009| >. Usually, efficient 
particle acceleration in these scenariosls related either to gap- 
type (jNeronov & Aharonian 2 007 1 [Levinson & Rieger|2011[ l 
or centrifugal-type processes ([Rieger & Aharonian 2008 ) oc- 
curring in the magnetosphere around a rotating black hole. 
While the former can be very efficient and lead to the onset of 
an electromagnetic pair cascade (triggered by the absorption 
of ambient photons up-scattered to high energies by electrons 


accelerating in the gap), the latter is less efficient such that the 
VHE spectrum is expected to be sh aped by the ambient soft 
photon spectrum (see also Vincent & LeBohec [2010 1 . Both 
approaches appear to be able to reproduce the observed hard 
VHE flaring characteristics, but may need an additional con- 
tribution to account for the HE continuum. 

Recently, a different type of modeling has been brought 
forward to explain the observed VHE y-ray emission from 
M87: Barkov et al. ( 2010) > argued that interactions of a rel- 
ativistic ancTstrongly magnetized outflow around the jet for- 
mation zone with a star partially tidally disrupted by the in- 
teraction with the SMBH can lead to a very efficient produc- 
tion of hadronic-ori ginating VHE photons (see also Bednarek| 


& Protheroe 1997 for a different version of the jet-star in- 
teraction model for blazar-type sources). The predicted ex- 
ponential character of the flux increase and decay during the 
resulting VHE event, a flat GeV to TeV y-ray spectrum from 
proton-proton interactions, as well as the accompanying X- 
ray flux enhancement due to the free-free emission of the 
shocked cloud of the stellar matter, are particularly interest- 
ing features of the model which should be kept in mind. 
Other scenarios could yet be considered in the same con- 
text as well, involving stochastic acceleration of high en- 
ergy particles within a turbulent accretion flow close to the 
event horizon of a SMBH (in analogy to the model proposed 
by Liu et al7]|2006] in the context of VHE y-ray emission 
of Sgr A*), or reconnection-driven impulsive acceleration of 
electron-positron pairs in a magnetized corona of the accre- 
tion disk, e.g., in analogy to the model proposed byjZdziarski 


|et al.|(|2009|) in the context of VHE observations of the Galac- 


tic X-ray binary system Cygnus X-l (see also de Gouveia Dal 
[Pino et al.|2010) . 


HST-1 — The arguments against association of the detected 
VHE y-ray emission with the non-thermal activity of the HST- 
1 knot are twofold: First, rapid variability of the VHE con- 
tinuum involving day-long flares implies that the emission 
region size is significantly smaller than the inferred size of 
the knot. Second, apart from the 2005 event, no correla- 
tion between the VHE activity and the radio-to-X-ray syn- 
chrotron radiation of the HST-1 flaring region, variable on the 
timescales of weeks and months, has been established. The 
first argument has to be taken with caution, though. That is 
because the HST-1 flaring region is basically unresolved even 
for radio interferometers, as already noted before. More im- 
portantly, very recently it has been found that rapidly variable 
high-energy radiation can be generated far away from the cen- 
tral jet engine, within the outer parts of relativistic outflows. 
In particular, the analysis of the y-ray emission of a bright 
quasar PKS 1222+216 indicates that the day-long VHE flares 
in this object originate most likely in compact sub-volumes 
of a relativistic jet at parsec distances from the active cen- 
ter (|Aleksic et al.||2011| |Tanaka et al.||201 f[ |Tavecchio et al.| 
20 1 1 [ . Interestingly, analogous phenomena seem to occur on 
even larger scales as well. For example, Chandra monitoring 
of the Pictor A radio galaxy reveals that the synchrotron X-ray 
emission of the knots located at tens of kpc from the jet base 
is variable on the timescales of years. This is much shorter 
than the variability timescale expected following the causality 
arguments for a given jet radius at the position of the knot, 
which is of the order of thousands of years (Marshall et al. 
2010|. Hence, it remains formally possible that also in the 
case of the HST-1 knot, located about 100 pc from the M87 
nucleus, rapid VHE flares are being generated on timescales 
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shorter than the ones expected for a homogeneous outflow. 

Yet, the indications found in this paper for correlation be- 
tween the VHE flares and an X-ray flux increases of the nu- 
cleus of M 87 provide observational support for the idea that 
the observed VHE y-ray emission, at least during the flaring 
episodes, is associated with the innermost parts of the jet or 
with the closest vicinities of the SMBH in the center of this 
radio galaxy. Again, in principal, it remains plausible that, 
while the unresolved core dominates the flaring states, the qui- 
escent VHE y-ray emission of M 87 is instead related to the 
HST-1 knot. However, above we have advocated for a com- 
mon origin of the quiescence and flaring VHE fluxes based 
on the spectral similarity of the two states. If the favored hy- 
pothesis is correct indeed, then interesting constraints on the 
physical parameters of HST-1 can be derived. That is because 
the inverse-Compton up-scattering of ambient photon fields 
to the VHE range by the high-energy electrons producing the 
observed radio-to-X-ray synchrotron radiation of the knot is 
at some level inevitable, with the efficiency depending most 
crucially on the unknown magnetic field intensity within the 
considered jet region. Since for the equipartition value of the 
jet magnetic field a relative ly intense VHE y-ray emission of 
HST- 1 should be expected (Stawarz et al.|2006 Cheung et al.' 


2007)>, the upper limits for the radiative output of the knot in 


y- rays (following from the fact that the observed VHE flux is 
associated with the nucleus) implies a strong, possibly even 
dynamically relevant jet magnetic field at hundreds-of-pc dis- 
tances from the jet base. This should be then considered as an 
important support for the MHD models of AGN jets, like the 
one presented in the particular context of the M 87 core and 
HST-1 knot by ( Nakamura et al.[2010) >, and discussed further 
from the observational perspective by |Perlman et al.| (|2011 1 
and|Chen et al.|(|201 1|). Analogous constraints emerging from 
the VHE data ca n be investigated for the outer parts of the 
M 87 jet as well (IStawarz et al. 1120051 IHardcastle & Croston 
|20TT) . 


7. SUMMARY & CONCLUSIONS 

During a joint VHE monitoring campaign on the nearby 
radio galaxy M87 in 2010, a major flux outburst was de- 
tected, triggering further VHE and MWL observations. The 
coordinated observations led to the best-sampled VHE light 
curve during a flaring state from this source (21 observa- 
tions in 15 days), revealing a single, isolated outburst. The 
measured VHE light curve of the flare is well described by 
a two-sided exponential function with significantly different 
flux doubling times of r™ e = ( 1 .69 + 0.30) days for the ris- 
ing and rj ecay = (0.611 + 0.080) days for the decaying part, 
i.e. a difference of a factor 2.77 + 0.62. This measurement 
provides the shortest and the most precisely determined VHE 
variability timescale of M 87 today, and the first detection of 
a significantly asymmetric VHE flare profile in the source. 

In comparison to previous VHE flares detected in 2005 and 
2008, the 2010 flare shows similar timescales and peak flux 
levels, but the overall variability pattern is somewhat different 
from the more extended periods of flaring activity with sev- 
eral flux maxima observed before, though the statistics and 
the sampling of the previous VHE flares limit a definitive con- 
clusion on that matter. From the VHE long-term light curve 
the duty cycle for VHE flares is estimated to be < 4 - 28%, 
depending on the assumed threshold flux defining a VHE high 
state. 

VLBA 43 GHz observations, triggered by the detection of 


the VHE flare, show no indications for an enhanced radio 
emission from the jet base in 2010. This is in contrast to ob- 
servations in 2008, where the detection of a radio outburst of 
the core contemporaneous with the VHE flare lead to the con- 
clusion that the VHE y-ray emission is likely produced in the 
direct vicinity of the SMBH (Acciari et al. 20091. Chandra 
X-ray observations, taken ~ 3 days after the peak of the VHE 
y-ray emission, show a high flux state of the core region in 
2010, supporting the interpretation that the VHE flare origi- 
nates from the innermost jet regions. 

The long term (2001-2010) light curve of M87, spanning 
from radio to VHE, is investigated for a common MWL signa- 
ture accompanying the three VHE flares. No unique signature 
is found. Observations of the jet component HST-1, which 
has also been proposed as a possible emission site of the VHE 
y-ray emission, show no enhanced activity in 2008 and 2010, 
disfavoring HST-1 as the origin of the VHE flares during those 
years. In 2008 and 2010 the VHE flares are accompanied by 
a high state of the X-ray core with a flux increase by a fac- 
tor ~ 2, while in 2005 the strong flux dominance (more than 
a factor 30) of the nearby X-ray feature HST-1 could have 
suppressed the detection of such an increase of the core emis- 
sion. Associating the VHE flares with the X-ray flares from 
the core places the emission site in the central resolution ele- 
ment of Chandra (0.6 " ~ 50 pc). 

Several models have been proposed to explain the observed 
VHE y-ray emission from M 87, most of which were based 
on the ‘misaligned BL Lac’ hypothesis. And, in fact, sev- 
eral observed properties of the source, including the broad- 
band y-ray spectrum, the multi-wavelength character of the 
VHE flares, and their overall energetics, could be consid- 
ered as support for the blazar-like models for the VHE y- 
ray emission. Yet, the particular one-zone or two-zone emis- 
sion models proposed involving only moderate beaming, both 
hadronic and leptonic, have difficulties in explaining VHE/X- 
ray correlated variability, for which indications are found in 
this paper, as well as the relatively flat VHE spectrum of 
M87 (photon index L ~ 2.3) extending up to ~ lOTeV. 
More complex blazar-type scenarios, involving ultra-compact 
emission regions and more substantial beaming, may, in- 
stead, face difficulties in accounting for the flares’ energet- 
ics with the isotropic peak luminosity during the flares reach- 
ing L V he — (0.8 - 2.4) x 1 0 42 erg s 1 , which is of order of 
the accretion luminosity and about 1% of the total jet power. 
Other non-blazar emission models discussed in the context 
of the VHE observations of M 87, which are quite success- 
ful in reproducing some of the established properties of the 
source, seem, on the other hand, challenged by the broad-band 
character of the y-ray continuum, as long as the GeV photon 
energy range constrained by the Fermi-LAI observations is 
included. All in all, even though no emission model could 
be identified as the most plausible one, important limitations 
for most of them have been identified and novel observational 
constraints have been presented in this paper. Thus, the case 
of M 87 demonstrates the relevance of systematic, long-term 
and multi-wavelength monitoring of nearby radio galaxies in 
understanding the origin of the high-energy emission of radio- 
loud AGN. 
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APPENDIX 

QUANTIFYING THE DIFFERENCE BETWEEN THE 2008 AND THE 2010 VHE FLARES 

To investigate the compatibility of the 2008 and 2005 VHE data with the flare detected in 2010, the best fit exponential function 
to the 2010 data (Fig. [3] Tab. [2} is taken as a template and compared with the 2008 and 2005 light curves. The position of the peak 
to is varied (all other parameters fixed) and the resulting ;y 2 /d.o.(. and probabilities are calculated for the data within a time span 
of length of the original fit length: ~ 7 < Gam < U) + 2. The results for the 2010, 2008, and 2005 data around the flaring episodes 

are shown in Fig.pl with the gray band again denoting the observation period with the flaring episode in each year. The statistics 
in 2005 are low resulting in several good probabilities (high probabilities > 0.1) over the flaring episode. For most parts of the 
2008 flare, the probability for the 2010 fit describing the data is well below 1%. The only notable exception is the beginning 
of the episode around MJD 54496 where the first two data points are reasonably well described by the function resulting in a 
probability of 0(5%). Note that the peak flux from the fit (~ 2 x 1CT 11 cirr 2 s _1 ) is similar to the peak fluxes measured in 2008. 
Varying the fit parameters in their errors does not significantly change the results. For comparison, the same study is shown for 
the 2010 data where, at the best fit position from the fit, the probability is ~70%. 
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